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Therapeutic  Mechanism  of  BET  Bromodomain  Inhibitor  in  Breast  Cancer 
Introduction 

Epigenetic  drugs  that  interact  with  chromatin  regulators  to  modify  the  chromatin 
states  of  cancer  cells  promise  to  be  important  in  cancer  treatment  (Yoo  and  Jones, 
2006;  Baylin  and  Jones,  2011).  To  effectively  use  these  drugs  to  manipulate 
transcription  for  therapeutic  purposes,  we  need  to  gain  insight  into  their  genome¬ 
wide  function.  We  propose  to  develop  both  the  experimental  and  computational 
methods  and  apply  them  to  identify  the  efficacy,  mechanism,  and  targets  of  BET 
domain  inhibitor  JQl  in  different  breast  cancer  cells.  This  could  potentially  inform 
personalized  breast  cancer  treatment  based  on  tumor  subtypes.  The  proposed 
experimental  and  computational  approaches  will  be  widely  applicable  to  the 
development  of  other  epigenetic  drugs  being  developed  to  treat  breast  cancers. 

Aim  1:  Identify  the  cistrome,  transcriptome  and  chromatin  dynamics  in  a  breast 
cancer  cell  line  SUM159  treated  with  BET  domain  inhibitor  JQl. 

Aim  2:  Computationally  infer  and  experimentally  validate  the  JQl  induced 
transcriptional  regulatory  network  and  pathway  in  SUM159. 

Aim  3:  Apply  the  above  approaches  to  test  JQl  on  different  breast  cancer  cell  lines 
and  xenografts  to  compare  the  efficacy  and  specificity  of  JQl  in  different  breast 
cancer  cells. 

Keywords 

Breast  cancer,  bromodomain  and  extra  terminal  domain  inhibitors,  JQl,  luminal 
breast  cancer,  basal  breast  cancer,  BRD4,  Chip-seq 

Overall  Project  Summary 

Task  1.  Identify  cistrome,  transcriptome  and  chromatin  dynamics  in  the  breast 
cancer  cell  line  SUM159  treated  with  BET  domain  inhibitor  JQl  or  BRD4 
knockdown. 

To  optimize  the  DNase-seq  protocol,  we  digested  2.5x10^  breast  cancer  MCF-7  cells 
with  50  units  of  DNase  1  (Roche),  and  then  selected  DNA  fragments  in  50-100,  100- 
200,  and  200-300bp  sizes  from  each  of  the  digested  samples  and  barcoded  to 
multiplex  HiSeq  sequencing.  We  carried  out  an  analysis  comparing  the  ability  of 
data  sets  from  different  size  fractions  to  discover  transcription  factor  binding  sites 
that  had  been  identified  by  ChlP-seq.  We  found  that  for  the  estrogen  receptor  and 
CTCF  short  fragments  were  optimal  (Fig.  1).  Our  analysis  of  the  location  of  DNasel 
cleavage  sites  with  respect  to  nucleosome  positions  revealed  that  chromatin 
structure  is  likely  to  be  the  reason  for  the  effectiveness  of  the  short  fragment  DNase- 
seq. 
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Figure  1.  Fraction  of  estrogen  receptor  (ER)  and  CTCF  ChIP-seq  peaks  in  MCF7 
identified  as  DNase-seq  peaks  for  different  fragment  size  selections.  Short  fragments 
are  much  more  efficient  than  longer  ones  in  recovering  CTCF  binding  sites  but  this 
effect  is  less  pronounced  for  the  estrogen  receptor. 


Preliminary  data  (Fig  2)  shows  that  DNase-seq  data  looks  highly  similar  to  BRD4 
ChlP-seq  and  did  not  provide  much  new  information  about  the  influence  of  JQl  on 
chromatin  and  transcriptional  changes.  This  is  a  consistent  with  the  known  role  of 
BRD4  in  opening  chromatin  (Voigt  and  Reinberg,  2011),  although  the  degree  of 
similarity  between  BRD4  ChlP-seq  and  DNase-seq  is  remarkable.  In  subsequent 
analyses  we  will  use  H3K27ac,  a  mark  that  is  indicative  of  enhancer  activity,  in 
addition  to  BRD4  ChlP-seq. 
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Figure  2.  BRD4  ChlP-seq  and  DNase-seq  in  untreated  MCF-7  cells  and  after  Ih  and 
6h  of  JQl  treatments.  BRD4  ChlP-seq  signal  is  very  similar  to  DNasel 
hypersensitivity  and  does  not  provide  much  complementary  information. 
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Task  2:  Computationally  infer  and  experimentally  validate  the  JQl  induced 
transcriptional  regulatory  network  and  pathway  in  SUM159. 

We  conducted  growth  rate  analyses  on  a  panel  of  breast  cancer  cell  lines  treated  by 
different  concentrations  of  JQl.  This  analysis  revealed  that  basal  breast  cancer  cell 
lines  tend  to  be  more  responsive  to  JQl  treatment  than  luminal  cell  lines.  For  this 
reason  we  decided  to  further  study  the  treatment  of  the  basal  SUM159  cell  line 
instead  of  the  luminal  MCF-7  cell  line. 

We  analyzed  the  genome-gene  expression  profile  of  SUM159  after  JQl  treatment  at 
3h,  12h  and  24h  in  relation  to  a  control  group  treated  with  DMSO.  This  analysis 
revealed  that  genes  were  predominantly  down-regulated  after  short  term  (3h) 
treatment  by  JQl,  as  expected  from  current  models  of  JQl  action  (Filippakopoulos  et 
al.,  2010).  We  used  a  model  of  "regulatory  potential"  to  use  BRD4  to  predict  the 
genes  likely  to  be  regulated  by  BRD4  and  effected  by  JQl  treatment.  In  this  model  a 
score  is  calculated  for  each  gene  as  a  sum  of  weights  in  which  BRD4  ChlP-seq  peaks 
close  to  the  transcription  start  site  of  the  gene  are  given  more  weight  than  BRD4 
peaks  that  are  far  from  the  TSS  of  that  gene.  We  evaluated  the  performance  of  this 
model  in  predicting  gene  expression  changes  associated  with  3h,12h  and  24h  of  JQl 
treatment  and  found  that  BRD4  ChlP-seq  could  inform  which  genes  were  down- 
regulated  at  3h  but  could  not  predict  the  up-regulated  genes.  The  performance  of 
this  model  deteriorated  with  treatment  time  as  secondary  effects  dominated  the 
primary  effect  of  BET  inhibition.  A  similar  analysis  based  on  H3K27ac  ChlP-seq  data 
resulted  in  a  poorer  prediction  for  down-regulated  genes  but  also  some  predictive 
power  for  up-regulated  genes. 


Figure  3.  Regulatory  potential  model  performance  using  BRD4  ChlP-seq  to  predict 
differentially  expressed  genes  for  3h  (A),  12h  (B),  and  24h  (C)  of  JQl  treatment. 
BRD4  ChlP-seq  can  predict  early  (3h)  down-regulated  genes  but  not  up  genes  or  late 
(24h)  genes. 

Task  3:  Apply  the  above  approaches  to  test  JQl  on  different  breast  cancer  cell  lines 
and  xenografts  to  compare  the  efficacy  and  specificity  of  JQl  in  different  breast 
cancer  cells. 
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We  have  collaborated  with  Dr  Polyak’s  lab  to  assess  the  efficacy  of  JQl  on  different 
breast  cancer  cells  lines  and  discovered  that  JQl  tends  to  slow  the  growth  of  basal 
cell  types  more  than  luminal  ones.  In  this  collaboration  we  have  also  assessed  the 
effect  of  long-term  treatment  of  SUM159  cells  with  JQl.  After  several  months  of 
treatment  we  find  these  JQl  sensitive  cells  become  less  sensitive  to  JQl.  RNA-seq 
analysis  of  this  JQl  insensitive  cell  line  reveals  that  the  gene  expression  pattern 
appears  to  be  more  luminal  than  the  parental  SUM  159  cells. 

Key  Research  Accomplishments 

•  Optimized  DNase-seq  conditions  and  developed  a  model  that  explains 
fragment  size  effects  in  DNase-seq  data. 

•  Analyzed  and  integrated  BRD4  ChlP-chip  and  expression  data  in  SUM159 
cell  line,  before  and  after  JQl  treatment. 

•  Conducted  H3K27ac  histone  mark  ChlP-seq  experiments  in  SUM159  cell 
lines. 

•  Established  a  JQl  insensitive  cell  line  from  a  parental  cell  line  that  is  JQl 
sensitive. 

Conclusion 

To  understand  the  mechanism  of  JQl  resistance  and  its  relationship  with  cell  type 
we  will  conduct  functional  genomics  and  further  epigenetic  analyses.  Specifically, 
we  will  perform  a  genome-wise  CR1SPR-Cas9  screen  of  genes  in  parental  and 
resistant  SUM159  cell  lines,  assessing  which  genes  support  or  hinder  cell 
proliferation  in  the  presence  and  absence  of  JQl.  We  will  compare  gene  expression 
changes  as  well  as  BRD4  and  H3K27ac  ChlP-seq  changes  between  parental  and 
resistant  cell  lines  with  the  genes  that  are  revealed  by  this  screen.  Based  on  the 
results  of  the  screen  we  will  assess  the  effect  of  gene  knock-outs  on  genome-wide 
gene  expression  to  examine  the  extent  to  which  these  knock-outs  can  modulated  the 
basal  to  luminal  transition  and  thereby  yield  insight  into  the  higher  sensitivity  of 
basal  breast  cancer  cell  types  to  the  BET  inhibitor  JQl. 
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